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Effects of error sources on ORION baseline vector determination are studied. Partial 
derivatives of delay observations with respect to each error source are formulated. Covari- 
ance analysis is then performed to estimate the contribution of each error source to base- 
line vector error. System design parameters such as antenna sizes, system temperatures 
and provision for dual-frequency (S- and X-band) operation are discussed. 


I. Introduction 

System design of an Operational Radio Interferometry 
Observing Network (ORION) is underway. The basic net- 
work will consist of a base station having a relatively large 
antenna and a mobile station having a smaller antenna. Nomi- 
nal operating frequencies will be at X-band with a span band- 
width of up to 400 MHz. Among the cost-effective parameters 
to be designed are the antenna sizes, system temperatures and 
the provision for S- and X-band operation for possible need of 
charged media calibration. To understand how these param- 
eters will affect the ultimate baseline measurement accuracy, 
the effects not only of these parameters but of all other 
possible error sources need to be studied. This article provides 
a covariance analysis to give an estimate of the effects of each 
individual error source. The ultimate baseline determination 
accuracy is estimated by taking the root-sum-square (rss) of 
the effects from all error sources. Baseline lengths of from 300 
km to 5000 km are to be considered. Both 3-D baseline vector 
accuracy and baseline length accuracy are to be estimated. 

II. Delay Observation 

The information for the solution of baseline vector is 
derived from delay observations, over the baseline, of extra- 
galactic radio sources (EGRS). The error-free information con- 
tent of a single observation can be written as 

r = -B cos 5 sin //A - i? cos 5 cos + /? sin 6 (1) 

x y Z ' ' 


where 5 = EGRS declination angle, HA = EGRS hour -angle at 
baseline center, B x , B y , B z = baseline components with B z 
parallel to the Earth’s spin axis, B x perpendicular to the meri- 
dian plane at baseline center and B y completing a right-hand 
Cartesian coordinate system as shown in Fig. 1 , and r is 
delay over the baseline in km. 

To determine the three baseline components at least three 
observations with non-coplanar lines-of-sight are needed, with 
orthogonal observations (perpendicular lines-of-sight) being 
optimum (in the sense of decoupling the three baseline compo- 
nents from one another). 

In general, observations are sensitive to such systematic 
parameters as epoch offset, frequency drift, etc. Such sensi- 
tivity can be alleviated by including these parameters as addi- 
tional estimated parameters (Ref. 1). Furthermore, it is well 
known that the sensitivity of a solution vector to observation 
error decreases as the number of independent observations 
increases. On the other hand, increasing the number of inde- 
pendent observations increases the total observation time and 
the data volume to be processed. Hence a compromise needs 
to be made. 

When the number of observations is greater than three, the 
optimum combination is such that the direction vectors of the 
EGRS being observed have equal spreads in all three orthogonal 
directions. On the other hand, the availability of EGRS within 
a viewing window is limited. Specifically, the 14 sources listed 
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in Table 1 have been found to have correlated flux densities 
higher than 4 Janskys at X-band over a 300-km baseline. 
Hence an optimum observation sequence is not always attain- 
able. A near optimum sequence, however, may be approached 
by proper scheduling. Table 2 is an observation sequence out 
of the 14 sources over a 24-hour period. Any 6-hour subset of 
this sequence provides a reasonable coverage of source position 
vectors in any three orthogonal directions. The sequence is 
in terms of the mean right ascension (RA) of the baseline 
center. The nominal observation time is 12 minutes except for 
the two sources having flux densities higher than 20 Janskys 
for which 4-minute observation time is allotted. A l°/sec 
antenna slew rate has been assumed without considering cable 
wrapup allowances. 
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where 0 is the latitude of the baseline center and £2 is the 
lateral orientation angle of the baseline from due east. Hence 
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III. Covariance Analysis 

The sensitivities of a delay observation to a perturbation of 
the baseline vector components are the partial derivatives of 
Eq. (1) with respect to these components: 
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It is felt that perturbations in the length L and in the lateral 
LT and vertical VT directions of a baseline would provide 
better physical meaning than those in B x , B y and B z . The 
partial derivatives of a delay observation with respect to these 
three physically meaningful components of a baseline can be 
derived from Eq. (2) with the aid of the following transforma- 
tion matrix, which is straightforward in view of Fig. 2: 
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-cos 6 (sin HA cos £2 + cos HA sin £2 sin 0) 
+ sin 5 sin f2 cos 0 
cos 5 (sin HA sin £2 - cos HA cos O sin 0) 
+ sin 5 cos £2 cos 0 
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As mentioned in Section II, effects of clock parameters can be 
alleviated by estimating these parameters. These include the 
epoch offset T q and as many as 12 segments of clock rate 
offsets, f v The sensitivities of a delay observation 

to these parameters can easily be formulated: Let f n applies to 
observations within the time segment / < t < t n \ then the 

sensitivities of a delay observation to these additional param- 
eters within tills time segment are 



and 

~ t t _ j) , i<n 

c(t~t n _j) , / = n (6) 

0 , i>n 

where c is the velocity of light. 

For the solutions of N estimated parameters from M 
independent delay observations, an N g X N information 
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matrix J - A T WA can be constructed the inverse of which is 
the computed covariance matrix 

p x = r l = {a t way 1 (7) 

Here, A is an M X N g sensitivity matrix with elements 
calculated according to Eqs. (4) - (6) for the M observations; 
W is an MXM diagonal weighting matrix with W mm = 
o~ 2 , the inverse of the error variance of the m ih delay 

observation. The diagonal elements of the computed covari- 
ance matrix (7) are the variances of the estimated parameters 
due to random errors in delay observations. Such errors 
include the effects of system noise and random tropospheric 
and ionospheric delay errors. Effects of systematic errors, such 
as atmospheric delay model errors, clock instability and source 
position errors, can be accounted for by including these 
parameters as “consider” parameters. The covariance matrix 
due to such systematic errors is given (Ref. 2) by 

P' = P A T WCP C t WAP (8) 

XX C X v ' 

where P is the covariance matrix of the N consider 

C C 

parameters and C is the sensitivity matrix to these parameters 
of the M delay observations. 

Platform parameter errors (polar motion, UT1) and antenna 
location error are directly related to baseline vector error 
independent of observation sequence. Hence no covariance 
analysis of such errors is necessary. 

Note that the effects of each individual error source on 
baseline vector components are of interest. These are calcu- 
lated one at a time by setting all other error sources to be zero. 


IV. Error Sources 

This section calculates the variances a 2 of a delay 
observation due to random error sources and the sensitivities 
of the observation to systematic error sources. The variances 
o 2 will be used to construct the weighting matrix W in 
Eq. (7) and the sensitivities will be used in the construction of 
sensitivity matrix C in Eq. (8). Covariance matrices can then 
be calculated. The errors in the three baseline components are, 
respectively, given by the square roots of the first three 
diagonal elements of the corresponding covariance matrix. The 
following paragraphs study the effects of each error source. 


A. System Noise 

Delay observation error due to the system noise in the 
receiver front ends at the two sites of a baseline can be derived 
from Ref. 3 as 1 


° T 2r rf 
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c = speed of light 
f span = span bandwidth, Hz 
D. = antenna diameter at site i, m 
T . = system noise temperature at site i, K 
B = channel bandwidth, Hz 
e. = antenna efficiency at site i 
S = correlated flux density of EGRS, Janskys 
T = observation (integration) time, sec 

This is a random error and its effects on baseline vector 
components are calculated from the computed covariance 
matrix of Eq. (7). 

B. Ionosphere 

Ionospheric error can be divided into two types, systematic 
and random. Systematic error is a mismodeling error which 
can be expressed as 

° T = e peak ( 10 ) 

where f(X.) and g(y.) are the solar-zenith-angle factor and 
elevation-angle factor of site /, and e peak is the modeling error 
of the ionospheric peak at zenith. The solar-zenith-angle factor 
can be represented by the empirical function 

f(X) = 0.2 + 0.8 cos 2 / 3 X (11) 


and the elevation-angle factor is given by (Ref. 4) 


Si-r) = m e + -R] cos 2 T ] l ' 2 - [(R e+ h x ) 2 

-R 2 cos 2 7] }l(h 2 - ftj) (12) 

1 A “fudge factor” of 0.8 has been included to account for possible 
imperfect bandpass overlap and data processing loss. 
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where R e = 6370 km, h l =215 km and h 2 = 454 km. 

Random ionospheric error is the result of fluctuations from 
the nominal values which can be represented by 

a T = \ \K XJfaJ +f(X 2 )g(y 2 )] e ' peak (13) 


where (A ///) 2 is the Allan variance of the clock at site i and 
/( N m ) is a discrete function, at sampling time At, simulating a 
flicker noise of unity Allan variance given by the following 
recurrence formula: 

f(j)= 1.95 /(/-I)- 0.95/0 -2) 


where e' peak is the standard deviation of the peak zenith 
fluctuation difference between the two sites. 

C. Troposphere 

Both dry and wet components of tropospheric error can be 
divided into systematic and random parts. Systematic error 
can be written as 
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where e z . is the zenith tropospheric calibration error at site i. 


with g. being the / th sample of a Gaussian random number of 
unity standard deviation and g. = 0 for /<1. For a 60-sec 
sampling time N becomes the time of observation in minutes 
rounded to the nearest integer. 


Random tropospheric error due to inhomogeneity can be 
written as 


E. Source Positions 


1 1 
+ 
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where e’ z is the 1-a standard deviation of the zenith 
on each site. 




(15) 

fluctuation 


The angular position error of any single EGRS can be 
considered to be fixed. The right ascension and declination 
components of this error are assumed to be uncorrelated to 
and independent of each other and of any other sources. 
Hence the effects of each component of a source should be 
treated as an independent consider parameter and yet the 
effects of all sources can be included in a single estimation. 


D. Clock Instability 

The effect of imperfect clocks on a delay observation is 
directly related to the difference between the drifts of the two 
clocks. It consists of an epoch offset and a relative drift. 
Epoch offset is an estimated parameter which is well isolated 
from baseline vector components and needs no further 
discussion. Clock drift effect is of a flicker random noise (as 
integrated over an observation time of ~12 minutes). It differs 
from a white random noise in that its Allan variance remains 
constant over all integration time instead of a linear decrease 
with integration time as for white random noise. Hence its 
effects cannot be treated as random effects. Instead, it is to be 
considered as systematic. The sensitivities of delay observa- 
tions to such error can be simulated by a simple flicker noise 
model: 



(16) 


The sensitivities of a delay observation to the right- 
ascension and declination errors of an EGRS can be formu- 
lated from Eq. (1): 


-r— - B cos 5 cos HA - B cos 5 sin HA 
9a * y 


(18) 

07 

gj- = B x sin 6 sin HA + B y sin 5 cos HA + B z cos 8 

F. Polar Motion and UT1 

The effects of polar motion error and of UT1 error on a 
baseline are independent of observation sequence and need no 
covariance analysis for their determination. It has been readily 
shown (Ref. 5) that x- and y -components of polar motion 
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error can be translated into the following baseline component 
errors: 


-B 

Z 
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As before, these can be translated into errors in baseline length 
and in lateral and vertical directions: 


= 0 


where a is the Earth’s radius at the pole. 
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These can be translated, with the aid of Fig. 2, into errors in - oj g B cos f2 cos 0 e UT1 

baseline length and in lateral and vertical directions: VT 


= 0 


Therefore, UTl error will affect only baseline orientation but 
not baseline length. 
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G. Antenna Location 

Antenna location may vary from one visit to another. Such 
error in repeatability translates directly into baseline error. 


H. Phase Accountability 


Therefore, the effects of polar motion error will be on baseline 
orientation only; baseline length will not be affected. 


This is a random phase variation in the receiver bandpass. 
Its effects on baseline are estimated in the same way as system 
noise effects. 


The effects of UTl error on baseline vector can be derived 
from Eq. (1). Let B x =B £ cos X and B y =B g sin X; then (1) 
can be rewritten as 

t = - B g cos 5 sin {HA +X) + B z sin S . ( 21 ) 

A perturbation in UTl of ey TI corresponds to a perturbation 
in HA of (jj g e UT] where co, is the Earth’s spinning speed. 
Such error will result in an error in X of equal magnitude but 
of opposite sign. In turn, an error in X of e x = - cj e e UT1 will 
result in errors in B and B of 

x y 


V. Results 

In this section, covariance analysis will be performed for 
the determination of a 300-km baseline at mid-latitude 
(0 = 35°) with a sequence of 32 delay observations. Such 
32-observation (~6-hr) sequence is picked out of Table 2. To 
examine the dependence on the observation sequence, the first 
of the 32 consecutive observations runs through Table 2. 
Baseline error due to each error source will then be shown as a 
function of the starting RA of the baseline center. The 
following assumptions are to be made: 

Baseline: 300 km at 0 = 35° 

Antenna diameters: 26 m (base) 

3 m (mobile) 


bB 
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Antenna efficiencies: 50% 

System temperatures: 40 K (base) 

130 K (mobile) 
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Channel bandwidth: 2 MHz 

Span bandwidth: 400 MHz 

EGRS flux densities: As in Table 1 

Observation times: As in Table 2 

Phase accountability: 0.2 cm each site (phase calibration 
assumed) 

Clock instability: A///= 1 0~ 1 4 (hydrogen maser) 

Zenith dry tropospheric error: 0.5 cm each site (systematic) 

Zenith wet tropospheric error 2 : O.S/v^Tcm each site 

(systematic) 

0.5/ y/2 cm each site 
(random) 

Zenith ionospheric error: 10 cm peak (systematic) 

2 cm difference per 1000 km 
(random) 

Polar motion error: 30 cm 

UT1 error: 0.6 msec 

EGRS angular position error: 0.01"/ sjl in RA and dec 

Antenna location error: 0.5 cm (mobile site only) 

Figures 3—11 display the results of covariance analysis. 
Here only the resultants of 3-D baseline errors are shown. 
Their mean values are summarized in Fig. 12 together with the 
effects due to platform parameter errors and antenna location 
error. Also shown in Fig. 12 are the corresponding errors in 
baseline length. It is observed that baseline length error 
accounts for only 30% of the total (3-D) baseline error. When 
baselines of different length are considered, it is found that 
errors due to ionosphere, source positions, polar motion and 
UT1 increase nearly linearly with baseline length. Such errors 
are drawn cross-hatched in Fig. 12. Errors due to other error 
sources essentially remain the same for all baseline lengths. 
(Tropospheric errors may increase slightly due to lower 
elevation angles for longer baselines.) 

VI. Discussions 

The above analysis provides an estimate of the effects of 
each possible error source on baseline vector determination. 
However, such information may not be of direct benefit to 
system design. In particular, one would need to look into the 
requirements on front-end system temperatures and antenna 
sizes to meet a given overall precision. These two parameters 
are among a few that are cost effective and controllable in the 


2 

Assuming water vapor radiometers are used. Bias errors will appear as 
clock epoch error which is to be estimated. 


ORION System design. Other error sources are limited by 
state-of-the-art capabilities: for instance, clock stability, media 
and platform parameter calibrations and the knowledge of 
EGRS angular positions. Both antenna sizes and system 
temperatures affect baseline vector determination through 
Eq. (9), the equation for system noise effect. The variation of 
such effect with antenna sizes and system temperatures will be 
studied, keeping all other effects unchanged. To account for 
reduced correlated flux densities over longer baselines and for 
full capability of 14 pairs of 2-MHz channels of the Mark III 
VLBI System, the following changes of parameters are 
assumed: 

EGRS correlated flux densities: 2 Janskys for all sources. 

Total channel bandwidth: 14X2 MHz per sideband. 

Observation times: 10 min. 

Figure 13 displays the effects of system noise as functions 
of antenna sizes. Three combinations of system temperatures 
are considered: 130K/130K, 40K/130K, and 40K/40K. 
The higher temperature refers to an uncooled receiver front 
end and the lower temperature to a cooled front end. The tick 
marks along the abscissa indicate typical combinations of 
antenna diameters in meters. The absolute maximum of the 
system noise effect is limited by the minimum fringe phase 
signal-to-noise ratio (SNR) below which phase tracking may 
fail. This corresponds to a baseline vector error of ~3.1 cm 
and a baseline length error of ~1.2 cm. This threshold is 
shown as a dashed line in Fig. 13. 

For a base station with D = 26 meters and T s = 40 K and a 
mobile station with D = 3 meters and T s = 130 K, the system 
noise effects are found to be ~0.8 cm on baseline vector and 
~0.3 cm on baseline length. These are comparable with the 
results shown in Fig. 12 for a 300-km baseline. It is observed 
that for such, and shorter, baseline lengths the effects of 
system noise are comparable with other error sources and 
consideration of antenna sizes and system temperatures may 
play an important role in overall baseline determination 
accuracy. For baselines longer than 300 km, those errors 
which are cross-hatched in Fig. 12 are to be scaled directly by 
baseline length 3 . Hence, for baselines 1000 km or longer, 
effects from ionosphere, source positions, polar motion and 
UT1 will dominate. It becomes relatively immaterial as to 
what combination of antenna sizes and system temperatures 
are used. This is illustrated in Figs. 14 and 15 for baseline 
vector error and baseline length error respectively. For a 
5000-km baseline, the overall errors are effectively constant 
regardless of antenna sizes and system temperatures. Unless 
the knowledge of source positions, polar motion and UT1 can 
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Error from troposphere will also increase slightly due to lower 
elevation angles for longer baselines. 
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be improved by a factor of ~3, consideration of antenna sizes 
and system temperatures is insignificant for baselines longer 
than 2000 km, as long as the signals are sufficiently strong to 
maintain a fringe phase SNR exceeding the minimum thresh- 
old. 


The above study has been based on observations over an 
integration time of 600 sec. Of course, tradeoffs can be made 
between antenna sizes and observation time and between 
system temperatures and observation time: A decrease in 
observation time by a factor of 2 corresponds to a decrease in 
the product of antenna diameters by a factor of \/2 or an 
increase in one of the system temperatures by a factor of 2. 
Therefore, antenna sizes greater than or system temperatures 
smaller than are needed to attain a specified precision allow a 
shorter integration time (and hence smaller data volume) per 
observation. Figures 13—15 can be used for integration times 


T other than 600 sec by simply moving the scale on the 
abscissa to the left by a factor (600 jT) 1 ! 2 

In regard to Fig. 12, the S/X calibration of ionospheric 
error needs to be mentioned. With the projected capabilities of 
source-position, polar-motion and UT1 calibrations assumed, 
the effects on baseline vector determination of these errors are 
much higher than that of ionospheric error (with Faraday 
rotation calibration). The dominance of these errors over 
ionospheric error will remain for all baseline lengths since they 
all increase with baseline length. Hence the need of S/X 
calibration does not seem justified unless the knowledge of 
source positions, polar motion and UT1 can be improved. This 
statement is not true in the case of baseline length determina- 
tion since polar motion and UT1 do not affect its accuracy as 
shown in the lower part of Fig. 12. Ionospheric error without 
S/X calibration becomes a major error souce in the determina- 
tion of baseline length 1000 km or longer. 
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Table 1. X-band VLBI extragalactic radio sources 


Correlated 
flux density 
over 300-km baseline, 3 


Source 

Right-ascension 

Declination 

Janskys 

3C 84 

49.12° 

41.33° 

44.6 

DA 193 

88.01 

39.81 

5.2 

P 0605-08 

91.40 

-8.57 

4.5 

OJ 287 

132.99 

20.30 

4.7 

4C 39.25 

140.98 

39.26 

7.8 

3C 273 

186.64 

2.33 

32.5 

3C 279 

193.40 

-5.52 

12.2 

3C 345 

250.32 

39.90 

6.8 

NR AO 530 

262.56 

-13.05 

5.1 

P 1741-038 

265.34 

-3.81 

5.6 

OV -236 

290.87 

-29.28 

6.8 

3C418 

309.28 

51.14 

5.6 

P 2134+004 

323.72 

.47 

8.0 

3C 454.3 

342.87 

15.88 

6.3 


Correlated flux density may vary with time; for the current analysis, 
however, these crude estimates can be used without introducing 
unrealistic error. 



Table 2. ORION 24-hour observation sequence 


Right-ascension 
of baseline center 


Source 

^obs 

hr 

min 

sec 

DA 193 

12 

0 

0 

0 

3C 84 

4 

0 

8 

29 

3C 454.3 

12 

0 

17 

31 

P 2134+004 

12 

0 

29 

55 

3C418 

12 

0 

42 

47 

3C 84 

4 

0 

51 

51 

DA 193 

12 

1 

0 

20 

P 2134+004 

12 

1 

14 

15 

3C 454.3 

12 

1 

26 

39 

3C418 

12 

1 

39 

23 

3C 84 

4 

1 

48 

27 

DA 193 

12 

1 

56 

56 

P 0605-008 

12 

2 

9 

44 

3C 84 

4 

2 

18 

47 

3C 454.3 

12 

2 

27 

49 

P 0605-08 

12 

2 

41 

39 

3C418 

12 

2 

55 

46 

3C 84 

4 

3 

4 

50 

P 0605-08 

12 

3 

13 

53 

4C 39.25 

12 

3 

26 

59 

3C 84 

4 

3 

36 

5 

DA 193 

12 

3 

44 

34 

OJ 287 

12 

3 

57 

17 

P 0605-08 

12 

4 

10 

7 

3C84 

4 

4 

19 

10 

DA 193 

12 

4 

27 

39 

OJ 287 

12 

4 

40 

22 

4C 39.25 

12 

4 

52 

42 

P 0605-08 

12 

5 

5 

48 

3C 84 

4 

5 

14 

51 

P 0605-08 

12 

5 

23 

54 

DA 193 

12 

5 

36 

42 

OJ 287 

12 

5 

49 

25 

P 0605-08 

12 

6 

2 

15 

3C 84 

4 

6 

11 

18 

DA 193 

12 

6 

19 

47 

P 0605-08 

12 

6 

32 

35 

4C 39.25 

12 

6 

45 

41 

OJ 287 

12 

6 

55 

1 

DA 193 

12 

7 

10 

44 

P 0605-08 

12 

7 

23 

32 

3C 84 

4 

7 

32 

35 

OJ 287 

12 

7 

41 

47 

4C 39.25 

12 

7 

54 

7 

P 0605-08 

12 

8 

7 

13 

3C 273 

4 

8 

16 

49 

3C 84 

4 

8 

22 

51 

OJ 287 

12 

8 

32 

3 

3C 84 

4 

8 

41 

15 

3C279 

12 

8 

51 

27 

3C273 

4 

8 

59 

37 

P 0605-08 

12 

9 

9 

13 

DA 193 

12 

9 

22 

1 

OJ 287 

12 

9 

34 

44 

4C 39.25 

12 

9 

47 

4 

3C273 

4 

9 

55 

59 

3C 279 

12 

10 

4 

9 

4C 39.25 

12 

10 

17 

15 

3C345 

12 

10 

30 

33 

3C 273 

4 

10 

39 

42 

OJ 287 

12 

10 

48 

37 

3C279 

12 

11 

1 

42 

DA 193 

12 

11 

15 

27 


Right-ascension 
of baseline center 


Source 

^obs 

hr 

min 

sec 

4C 39.25 

12 

11 

28 

7 

3C 273 

4 

11 

37 

2 

DA 193 

12 

11 

46 

37 

3C 345 

12 

12 

0 

16 

3C279 

12 

12 

13 

25 

3C 273 

4 

12 

21 

35 

OJ 287 

12 

12 

30 

30 

3C 279 
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Fig. 1. Geometry of baseline and EGRS 


Fig. 2 


B x = Bcc*ft dB L = dB (LENGTH) 

B = BsinA cos £ dB^. = B dTL (LATERAL) 

B^ = B sin sin £ dByj = B sin ild£ (VERTICAL) 

tt/2 - i = 4> = BASELINE CENTER LATITUDE 

. Transformation of baseline component differentials 























Fig. 11. Effects of system phase accountability error on baseline vector determination 
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Fig. 12. Baseline vector and length determination errors 
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Fig. 13. Contribution of system noise to baseline determination error 


30 




Fig. 14. Total baseline vector determination error 



Fig. 15. Total baseline length determination error 
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